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ABSTRACT 

Located in the Perseus cluster, NGC 1271 is an early-type galaxy with a small effective radius of 2.2 kpc 
and a large bulge stellar velocity dispersion of 276 km s“^ for its X-band luminosity of 8.9 x 10'° Lq. We 
present a mass measurement for the black hole in this compact, high-dispersion galaxy using observations 
from the integral field spectrograph NIFS on the Gemini North telescope assisted by laser guide star adap¬ 
tive optics, large-scale integral field unit observations with PPAK at the Calar Alto Observatory, and Hubble 
Space Telescope WFC3 imaging observations. We are able to map out the stellar kinematics both on small 
spatial scales, within the black hole sphere of influence, and on large scales that extend out to four times the 
galaxy’s effective radius. We find that the galaxy is rapidly rotating and exhibits a sharp rise in the velocity 
dispersion. Through the use of orbit-based stellar dynamical models, we determine that the black hole has a 
mass of (3.0^}'[) x 10^ Mq and the H-band stellar mass-to-light ratio is 1.40;|;q}j Tq (Ict uncertainties). NGC 
1271 occupies the sparsely-populated upper end of the black hole mass distribution, but is very different from 
the Brightest Cluster Galaxies (BCGs) and giant elliptical galaxies that are expected to host the most massive 
black holes. Interestingly, the black hole mass is an order of magnitude larger than expectations based on the 
galaxy’s bulge luminosity, but is consistent with the mass predicted using the galaxy’s bulge stellar velocity 
dispersion. More compact, high-dispersion galaxies need to be studied using high spatial resolution observa¬ 
tions to securely determine black hole masses, as there could be systematic differences in the black hole scaling 
relations between these types of galaxies and the BCGs/giant ellipticals, thereby implying different pathways 
for black hole and galaxy growth. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: individual (NGC 1271) - galaxies: kine¬ 
matics and dynamics - galaxies: nuclei - black hole physics 


1. INTRODUCTION 

Over the past 15 years, it has become increasingly clear 
that supermassive black holes are an essential component of a 
galaxy, as demonstrated by the correlations connecting black 
hole masses and galaxy bulge properties (e.g., Kormendy & 
Richstone 1995; Ferrarese & Merritt 2000; Gebhardt et al. 
2000; Marconi & Hunt 2003; Giiltekin et al. 2009; Kormendy 
& Ho 2013). Supermassive black holes are thought to regu¬ 
late galaxy properties and influence star formation via feed¬ 
back processes (Silk & Rees 1998; Fabian 1999), however 
the black hole - bulge relations can also arise because of the 
inherent averaging associated with random galaxy mergers, 
without the need for the black hole to actively influence its 
host galaxy (Peng 2007; Jahnke & Maccio 2011). Roughly 80 
dynamical black hole mass (Mbh) measurements have been 
made to date (Kormendy & Ho 2013), almost exclusively 
through the use of high angular resolution facilities such as 
the Hubble Space Telescope (HST) and 8-10m ground-based 
telescopes with the aid of adaptive optics (AO). Despite the 
growing number of black hole mass measurements, the local 
black hole mass census is highly incomplete. Gaining a more 
complete picture of black hole demographics and a deeper un¬ 
derstanding the mechanisms that drive black hole/galaxy evo¬ 
lution requires the secure measurement of many more black 
holes, particularly those at the extremes of the black hole mass 


scale and in a wider range of galaxy types. 

With the goal of finding more objects suitable for future 
dynamical black hole mass measurements, the HET Massive 
Galaxy Survey obtained long-slit spectra of ^1000 galaxies 
using the Hobby-Eberly Telescope (HET) at McDonald Ob¬ 
servatory (van den Bosch et al. 2015). The survey uncovered a 
sample of interesting early-type galaxies with small sizes and 
high stellar velocity dispersions for their luminosities. More 
quantitatively, these galaxies have an effective radius, Rg, be¬ 
low 3 kpc, a central stellar velocity dispersion larger than 
250 km s“*, and K-band luminosities ~(5-25) x 10^° Lq. 
The HET spectra hint that these compact, high-dispersion 
galaxies could host some of the largest black holes known 
(Mbh > 10® Mq), and that the black holes could weigh a 
high fraction of its host galaxy’s mass. Six example ob¬ 
jects were highlighted in van den Bosch et al. (2012), and 
orbit-based stellar dynamical models were calculated for one 
galaxy, NGC 1277. In the case of NGC 1277, van den Bosch 
et al. (2012) found a 1.7 x 10^® Mq black hole that is surpris¬ 
ingly 59% of the galaxy’s bulge mass, or 14% of the galaxy’s 
total mass. While Yildirim et al. (2015) infer a similar Mbh 
from seeing-limited, large-scale integral field unit (IFU) data, 
Emsellem (2013) show a smaller black hole of a few billion 
solar masses can also reasonably reproduce the observed HET 
long-slit kinematics presented in van den Bosch et al. (2012). 
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Obtaining secure black hole mass measurements for the 
other compact, high-dispersion galaxies found through the 
HET Massive Galaxy Survey is important for addressing 
questions concerning the upper end of the black hole - host 
galaxy relationships. With the present sample of black hole 
mass measurements, the slope, intrinsic scatter, and even the 
shape of the correlations for high-mass black holes are not 
well established (e.g., McConnell & Ma 2013). Also, the 
black hole mass - stellar velocity dispersion relation (Mbh - 
cr*) and the black hole mass - bulge luminosity relation (Mbh 

- Lbui) are in direct conflict at the upper end, and make drasti¬ 
cally different predictions for the inferred number density of 
the most massive black holes (Lauer et al. 2007). Not only 
are the compact, high-dispersion galaxies from the HET sur¬ 
vey useful for filling in the poorly sampled high-mass end of 
the black hole relations, but also the mass estimates from Mbh 

- CT* and Mbh - Tbui differ by a factor of at least three. Con¬ 
sequently, the galaxies are also useful for testing which of 
the correlations is more fundamental and a better predictor of 
Mbh at the high-mass end of the scaling relations. 

While recent progress has been made in searching for, and 
revising measurements for, black holes with masses larger 
than 10® Mq (Shen & Gebhardt 2010; Gebhardt et al. 2011; 
McConnell et al. 2011, 2012; van den Bosch et al. 2012; 
Walsh et al. 2010, 2013; Rush et al. 2013), many of these 
galaxies are giant ellipticals or Brightest Cluster Galaxies 
(BCGs), which are often large (with /?« > 10 kpc; e.g., Dalla 
Bonta et al. 2009), have cored surface brightness profiles, and 
are dispersion-supported showing little to no rotation. In con¬ 
trast, the compact, high-dispersion galaxies found through the 
HET survey are small, rapidly rotating, and generally exhibit 
cuspy surface brightness profiles. Such host galaxy environ¬ 
ments haven’t been extensively explored on the black hole - 
host galaxy relationships. Besides NGC 1277, only the com¬ 
pact galaxies NGC 1332 (Rush et al. 2011), NGC 4342 (Cret- 
ton &. van den Bosch 1999), NGC 4486B (Kormendy et al. 
1997), and M60-UCD1 (Seth et al. 2014) have dynamical 
black hole mass measurements, with NGC 1332 and NGC 
4342 being most like NGC 1277. Eor these galaxies, the 
black hole mass measurements are in agreement with Mbh 

- cr* given the intrinsic scatter of the relation, but are positive 
outliers from the Mbh - Tbui relation (Kormendy & Ho 2013). 
We note that there are uncertainties associated with the bulge 
luminosity for NGC 1332 (see Kormendy &. Ho 2013 for de¬ 
tails) and the black hole mass measurement for NGC 4486B 
(see Giiltekin et al. 2009 for details). Also, tidal stripping 
is believed to be the cause of the over-massive black hole in 
the ultracompact dwarf galaxy M60-UCD1 (Seth et al. 2014), 
and there is some debate as to whether NGC 4486B and NGC 
4342 have been stripped as well (e.g., Eaber 1973; Bogdan 
et al. 2012; Blom et al. 2014). Nevertheless, additional sim¬ 
ilar galaxies need to be studied because there could be sys¬ 
tematic differences in the scaling relations between the com¬ 
pact, high-dispersion galaxies and the giant ellipticals/BCGs. 
If so, that would imply that the black holes in the two types of 
galaxies grew in different ways. 

While the compact, high-dispersion galaxies are unusual in 
the present-day Universe, they are qualitatively similar to the 
typical z ~ 2 quiescent galaxies, which are also small, have 
disk-like features, and could have high velocity dispersions 
(Zirm et al. 2007; van Dokkum et al. 2009; van der Wei et al. 
2011). The z ~ 2 red nuggets are believed to be the progeni¬ 
tors of the massive early-type galaxies seen today, evolved in 


size and mass (e.g., van Dokkum et al. 2010). Thus, the com¬ 
pact, high-dispersion galaxies found through the HET survey 
could provide clues to the link between local galaxies and the 
z ^ 2 red nuggets, the ultra compact sub-mm galaxies at z ~ 3 
(Toft et al. 2014), and the early massive black holes found in 
z > 6 quasars (Eabian et al. 2013). 

We have begun to obtain the imaging and spectroscopic ob¬ 
servations necessary for a detailed examination of the com¬ 
pact, high-dispersion galaxies from the HET Massive Galaxy 
Survey. This includes HST and AO-assisted lEU observa¬ 
tions to probe the region over which the black hole domi¬ 
nates the galaxy’s potential (the black hole sphere of influ¬ 
ence; Tsphere = GMbh/ct*), and lEU observations that sample 
the large-scale stellar kinematics out to several effective radii. 
In this paper, we focus on measuring the mass of the black 
hole for the first compact, high-dispersion galaxy for which 
we have completed AO lEU observations. NGC 1271 has not 
been widely investigated in the literature and is given an un¬ 
certain SBO classification according to the NASA/IPAC Ex- 
tragalactic Database (NED). The galaxy is located within the 
Perseus cluster, at z = 0.0192, and we adopt a distance to NGC 
1271 of 80 Mpc, which is the Hubble flow distance derived 
from the Mould et al. (2000) Virgo + Great Attractor + Shap- 
ley Supercluster Infall velocity field model assuming a Hub¬ 
ble constant of Hq = 70.5 km s“' Mpc“', a matter density of 
VIm = 0.27 and a cosmological constant of Ha = 0.73. The 
Sloan Digital Sky Survey g-/ color is 1.6 and absolute r-band 
magnitude is -20.8 for the galaxy. Long-slit spectra of NGC 
1271 were obtained through the HET Massive Galaxy Survey, 
and the reported [N II]/Hq; and [O III]/H/3 emission-line ratios 
measured within a 3 "5 aperture (van den Bosch et al. 2015) 
place NGC 1271 just within the composite galaxies section 
near the active galactic nuclei side of Kewley et al. (2006). 
Also, van den Bosch et al. (2015) do not find the presence of 
any broad emission lines. 

In Section 2, we describe the imaging and spectroscopic 
observations, including the data reduction procedures. We 
present the luminous mass model for the galaxy in Section 
3 and the stellar kinematics in Section 4. Our determination 
of the point-spread function (PSE) for the spectroscopic ob¬ 
servations is discussed in Section 5, and an overview of the 
orbit-based stellar dynamical models is given in Section 6. In 
Section 7, we present the results of the models and examine 
possible sources of systematic uncertainty. Einally, in Sec¬ 
tions 8 and 9, we study the galaxy’s orbital structure, discuss 
NGC 1271 within the context of the Mbh - host galaxy rela¬ 
tionships, and summarize our findings. 

2. OBSERVATIONS 

Eor NGC 1271, we obtained imaging observations with the 
HST Wide-Eield Camera 3 (WEC3) in order to measure the 
galaxy’s surface brightness distribution. We also acquired 
spectra with the Near-infrared Integral Eield Spectrometer 
(NILS; McGregor et al. 2003) on the 8.1m Gemini North tele¬ 
scope assisted by the ALTtitude conjugate Adaptive optics for 
the InfraRed (Herriot et al. 2000; Boccas et al. 2006) system. 
The NILS data is important for constraining Mbh, as it re¬ 
solves the black hole sphere of influence. We also acquired 
large-scale spectra with the Postdam Multi Aperture Spec¬ 
trograph (PMAS; Roth et al. 2005) in the Pmas fiber PAcK 
(PPAK; Verheijen et al. 2004; Kelz et al. 2006) mode at the 
3.5m telescope at Calar Alto Observatory. Although long-slit 
spectroscopic observations along the galaxy major axis have 
been previously made using the HET, measuring a large-scale. 
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two-dimensional (2D) velocity field is preferable over a sin¬ 
gle slit observation for constraining the stellar mass-to-light 
ratio and the stellar orbital distribution. Hence, we use the 
PPAK IFU observations in place of the major-axis HET mea¬ 
surements. Below we describe the WFC3, NIFS, and PPAK 
observations and data reduction methods. 

2.1. HST Imaging 

We observed NGC 1271 with HST WFC3 and the 
IR/F160W filter under program GO-13050. The observation 
was composed of three dithered full array exposures of 450 
s, and four dithered subarray exposures of 1.7 s, leading to a 
total integration time of 1354 s. The short subarray exposures 
were chosen to ensure the nucleus would not become satu¬ 
rated. The flattened, calibrated images were corrected for geo¬ 
metric distortions, cleaned, and combined using AstroDrizzle 
(Gonzaga et al. 2012). Since the exposures were dominated 
by galaxy light, we found that the standard AstroDrizzle sky 
subtraction overestimated the background flux. We therefore 
manually measured the background level in each of the im¬ 
ages. For the full array images, we measured the flux from the 
corners of each image, while for the subarray exposures we 
measured the flux difference between the sky-subtracted full 
frames and the subarray frames. With the background level 
determined, the exposures were combined to produce a super- 
sampled image with a spatial resolution of 0."06 pixeF'. Not 
only is the HST image suitable for determining the luminous 
distribution on near the black hole, but due to the small size 
of NGC 1271, we were also able to measure the luminosity 
out to larger galaxy scales, extending to ~5 Re (adopting an 
effective radius of 5 "6, or 2.2 kpc, measured from a single 
component Sersic fit to the HST F160W image; see Section 
8 . 2 ). 

2.2. NIFS Spectroscopy 

The NIFS laser guide star (LGS) AO observations were ac¬ 
quired over three nights, on 2012 Dec 27, 2012 Dec 29, and 
2013 Jan 8, in queue mode under program GN-2012B-Q-51. 
We used the H+K filter and the K grating with a central wave¬ 
length of 2.2 pm to obtain spectra over a 3" x 3" field-of-view 
and a spectral resolution of R~5290. We recorded 900 s expo¬ 
sures of the galaxy nucleus, following an Object-Sky-Object 
observing sequence, totaling 3 hours of on-source integration. 
An R = 17.8 mag star located 17" away from the galaxy was 
used as the tip-tilt reference. In addition, we observed the tip- 
tilt star to monitor the PSF during each of the three nights and 
the AO V stars HIP 10559 and HIP 22842 for telluric correc¬ 
tion. 

The data were reduced using IRAF^ tasks within the Gem- 
ini/NIFS package version 1.11, utilizing the example NIFS 
processing scripts*. The reduction included sky subtraction, 
flat fielding, interpolation over bad pixels, cosmic-ray clean¬ 
ing, and spatial rectification and wavelength calibration using 
Ronchi mask and arc lamp exposures. The spectra were then 
corrected for telluric features, using an AO V star spectrum, 
after interpolating over the Br 7 absorption line and dividing 
by a black body with a temperature of 9480 K. Next, a data 
cube was produced, having x and y spatial dimensions, each 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron¬ 
omy under cooperative agreement with the National Science Foundation 

http://www.gemini.edu/sciops/instruments/nifs/data-format-and- 


with a scale of 0"05 pixel"', and one spectral dimension. A, 
using a common wavelength range and sampling for the in¬ 
dividual science exposures. The relative spatial positions be¬ 
tween the data cubes were determined by summing along the 
wavelength axis and cross-correlating the resulting flux maps. 
The offsets were used to align and combine the 12 individual 
exposures, generating the final data cube of the galaxy. Data 
reduction of the PSF star observations followed a similar pro¬ 
cedure. 

2.3. PPAK Spectroscopy 

The PPAK observations of NGC 1271 were acquired as part 
of a campaign to obtain large-scale spectroscopy of the com¬ 
pact, high-dispersion galaxies. The observations, data reduc¬ 
tion, and kinematic measurements for NGC 1271 will be pre¬ 
sented in Yildirim et al. (in prep), but follow closely the PPAK 
observations for two other compact, high-dispersion galax¬ 
ies described in Yildirim et al. (2015). For completeness, we 
briefly review the pertinent information below and in Section 
4.2. 

The wide-held IFU observations were taken over three 
nights, from 2013 Jan 4-6, using the V500 grating to pro¬ 
vide coverage of 4200-7000 A with a spectral resolution of 
R~850 at 5000 A. We used three dithers to fully sample the 
331 2."7-diameter hbers, to increase the spatial resolution of 
the data, and to address effects of vignetting and bad pixels. 
An additional 36 hbers, located 72" away from the center of 
the instrument held-of-view, were used to measure the sky. 
During each of the three nights, two 1200 s science exposures 
were taken at each of the three dither positions, leading to a 
total of 6 hours of on-source integration. The data reduction 
followed the procedure adopted for the Calar Alto Legacy In¬ 
tegral Field Spectroscopy Area Survey. The main steps in¬ 
cluded bias subtraction, Hat-helding, cosmic ray cleaning, ex¬ 
traction of spectra, wavelength calibration, sky subtraction, 
and Hux calibration using spectrophotometric standard stars. 
Spectra from the three pointings were then combined and re¬ 
sampled into a data cube, followed by a correction for dif¬ 
ferential atmospheric refraction. We note that the line spread 
function is measured as a function of wavelength and hber po¬ 
sition from arc lamps during the wavelength calibration step, 
and then homogenized to a common value prior to the extrac¬ 
tion of the PPAK kinematics. The details of these steps are 
discussed at length by Sanchez et al. (2012) and Husemann 
et al. (2013) and we refer the reader to those publications for 
additional information. 

3. CONSTRUCTING THE LUMINOUS MASS MODEL 

We generated a luminous mass model for NGC 1271 by pa¬ 
rameterizing the HST WFC3 F160W image as the sum of 2D 
Gaussians using the Multi-Gaussian Expansion (MGE) for¬ 
malism (Monnet et al. 1992; Emsellem et al. 1994). The MGE 
method is able to reproduce a wide range of galaxy surface 
brightness profiles and allows for an analytic deprojection to 
determine the intrinsic luminosity density. Here, we use the 
image decomposition package Galfit (Peng et al. 2010) be¬ 
cause it takes into consideration an error map during the fit 
and allows for the detailed examination of model residuals, 
but we utilize the implementation of Cappellari (2002) to de¬ 
termine suitable starting parameter values for the initial run 
with Galfit. When constructing the MGE model, we account 
for the WFC3 PSF, which we adopt from van der Wei et al. 
(2012). This PSF was generated with Tiny Tim (Krist &. Hook 
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Fig. 1.— Isophotes of the MGE model (red) are compared to the HST 
WFC3 F160W image (top) and to the inner 3”6x3”6 region of the image 
(bottom). Contours are logarithmically spaced, but arbitrary. Foreground 
stars and galaxies were masked during the MGE fit. 


2004) for the F160W filter at the center of the WFC3 detec¬ 
tor assuming a G2 V spectral type, and drizzled to produce 
a PSF with the same spatial scale as our final science image 
of NGC 1271. In addition, we identified foreground objects 
using the program SExtractor (Bertin & Arnouts 1996) and 
masked these regions during the MGE fit. 

The final MGE model contained 11 components, where 
each Gaussian was set to have the same position angle and 
center and the projected axis ratio {q') was required to be 
larger than 0.25. Due to the degeneracy associated with fitting 
a large number of Gaussians, we chose to restrict q' > 0.25 
in order to avoid highly flattened components that place very 
stringent constraints on the viewing angles for which a model 
can be deprojected. The final MGE model is an good repre¬ 
sentation of the galaxy, as can be seen in Eigure 1. We present 
the final MGE parameters, after correction for galactic extinc¬ 
tion using the Schlafly & Einkbeiner (2011) value given by 
NED, in Table 1. 

4. EXTRACTING THE STELLAR KINEMATICS 


TABLE 1 

MGE Parameters 


j 

(1) 

log Ij (Lq,// pc 2) 

(2) 

(t' (") 

(3) 

(4) 

I 

6.074 

0.073 

0.25 

2 

4.778 

0.175 

0.83 

3 

4.786 

0.399 

0.56 

4 

4.183 

0.765 

0.77 

5 

3.172 

1.838 

0.66 

6 

3.756 

2.018 

0.26 

7 

3.391 

4.306 

0.25 

8 

3.054 

5.857 

0.38 

9 

2.474 

8.876 

0.49 

10 

1.879 

12.946 

0.72 

11 

0.889 

24.772 

0.99 


Note. — The component number is listed in column (1), the cen¬ 
tral surface brightness, using a galactic extinction of 0.085 and a solar 
absolute magnitude of 3.33, is provided in column (2), the dispersion 
along the major axis is given in column (3), and the axis ratio is pre¬ 
sented in column (4). The components all have of a position angle of 
-50.3° and projected quantities are denoted with primed variables. 

Erom the NIES and PPAK data cubes, we measured the 
line-of-sight velocity distribution (LOSVD) as a function of 
spatial location. The LOSVD was described using the first 
four Gauss-Hermite (GH) moments: the radial velocity (V), 
the velocity dispersion (a), and /13 and /i 4 , which describe 
the LOSVD’s asymmetric and symmetric deviations from a 
Gaussian. High signal-to-noise (S/N) spectra, typically > 30, 
are required in order to reliably extract the higher order GH 
moments (e.g., van der Marel & Eranx 1993; Bender et al. 
1994), and thus we use the Voronoi binning algorithm (Cap- 
pellari & Copin 2003) in order to construct spatial bins that 
optimize the balance between spatial resolution and S/N. We 
then used the penalized pixel fitting (pPXE) method of Cap- 
pellari & Emsellem (2004) to measure the stellar kinemat¬ 
ics in each bin. This procedure determines the best-fitting 
LOSVD by convolving with a stellar template to match the 
observed galaxy spectra. Errors on the kinematics were de¬ 
termined using Monte Carlo simulations, in which random 
Gaussian noise was added to the spectrum based upon the 
pPXE model residuals. We performed 100 realizations and 
from the distributions measured the standard deviation to de¬ 
termine la uncertainties. During the Monte Carlo runs, the 
penalization term was set to zero to produce realistic errors. 

4.1. NIFS Kinematics 

We measured the stellar kinematics from the three pri¬ 
mary X-band CO bandheads [(2-0)'^CO, (2- l)'^CO, and 
(4-2)'^CO] in 127 spatial bins by using pPXE to fit the wave¬ 
length region between 2.26-2.42 ^m. We made use of the 
NIES Spectral Template Library v2.0® (Winge et al. 2009), 
which contains 28 stars observed using the NIES lEU with 
the K grating and H + K filter. The library includes spectral 
types ranging from G 8 - M5 giant stars, K3 - M3 supergiants, 
and a G 8 II star. 

We first created an optimal stellar template by fitting a high 
S/N spectrum, constructed by adding together all spectra in 
the galaxy data cube. This optimal template was composed 
of six stars, and was dominated by M5 III, M3 III, and K3 
lab stars that make up 37%, 26%, and 21% of the total flux, 
respectively. Next, we measured the GH moments in each 

® http://www.gemini.edu/sciops/instraments/nearir-resources/spectral- 
templates 
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spatial bin with pPXF by keeping the relative weights of the 
stars that make up the template fixed, but allowing the co¬ 
efficients of a second degree additive Legendre polynomial 
and a second degree multiplicative Legendre polynomial to 
vary. Such polynomials are needed to account for differences 
in the optimal stellar template and the galaxy spectra shape, 
as the continuum of the stars in the NIFS Spectral Template 
Library has been previously removed. The kinematics were 
in good agreement with those measured when fitting a new 
optimal stellar template to each spatial bin and fitting only 
two GH moments. Finally, we bi-symmetrize the kinematics 
using the machinery presented van den Bosch & de Zeeuw 
(2010). Since the dynamical models are only able to pro¬ 
duce symmetric kinematics, this step is commonly performed 
(e.g., Gebhardt et al. 2003; Cappellari et al. 2006; Onken et al. 
2014) in order to reduce the noise in the observations. With 
the symmetrization routine, the systematic offsets in the odd 
GH moments, such as the galaxy’s recession velocity, are re¬ 
moved as well. 

From the NIFS data, we find that the galaxy is rotating 
quickly, with stars reaching velocities of ±226 km s“'. There 
is also a sharp peak in the velocity dispersion, which rises 
from 205 km s“' at a radius of ^1" to 396 km s“' at the cen¬ 
ter. The map of /13 is anti-correlated with the velocity map, 
while hn shows a slight increase at the center. The S/N in each 
spatial bin (measured as the ratio between the median value of 
the spectrum and standard deviation of the pPXF model resid¬ 
uals) ranged between 33 and 96 with a median value of 66 . 
Therefore, we were able to place excellent constraints on the 
kinematics, with median errors over all spatial bins of 7 km 
s“', 9 km s“', 0.02, and 0.02 for V, a, h^, and hii, respectively. 
We present example spectra and fits with pPXF at three dif¬ 
ferent locations within the NIFS data cube in Figure 2 and the 
bi-symmetrized NIFS kinematics in Table 2. 

4.2. PPAK Kinematics 

We measured the stellar kinematics in 268 spatial bins over 
a wavelength range of 4200-7000 A, which includes a num¬ 
ber of absorption features such as the H/3, Mg Ib, and Fe 
5015 lines. During the fit with pPXF, we masked sky features 
and emission lines and included a 15th degree additive poly¬ 
nomial. The kinematics were extracted using the Indo-U.S. 
Library of Coude Feed Stellar Spectra (Valdes et al. 2004), 
and the optimal stellar template was dominated by G9 V, G9 
III, KO III, and AOp stars. As a final step, we bi-symmetrize 
the kinematics and subtract off the systematic offsets in the 
odd GH moments using the procedure in van den Bosch & de 
Zeeuw (2010). 

The kinematics from the PPAK data were measured out to 
^24", or ^4 Re- The large-scale kinematics exhibit features 
that are similar to the measurements made from the high spa¬ 
tial resolution NIFS data. In particular, the stars show rotation 
with velocities of ±231, a peak in the velocity dispersion to 
values of 297 km s“* from 102 km s“’ at radius of ^24", and 
there is a /13 - V anti-correlation. The difference in the peak 
velocity dispersions measured from the PPAK and NIFS data 
can be attributed to the very different spatial resolutions of the 
two data sets. Typical errors on the PPAK kinematics for V, cr, 
/ 13 , and are 8 km s“', 12 km s“', 0.04, and 0.05. We present 
example fits to the galaxy spectra at several spatial locations 
in Figure 3 and provide the bi-symmetrized PPAK kinematics 
in Table 3. 

NGC 1271 has an uncertain SBO classification according to 
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Fig. 2.— Shown in the top, middle, and bottom panels, respectively, are 
example NIFS spectra extracted from three spatial locations: a single spaxel 
located near the nucleus, a bin containing five spaxels at an intermediate dis¬ 
tance from the galaxy center, and one of the outermost bins composed of 28 
spaxels. Overplotted in red is the optimal stellar template convolved with the 
best-fitting LOSVD. The model residuals are shown in green, and have been 
shifted by an arbitrary amount. 

NED, but we do not see an obvious bar feature in the HST 
F160W image or in the PPAK/NIFS kinematics. N-body sim¬ 
ulations have shown that common kinematic signatures asso¬ 
ciated with bars include a “double-hump” feature in the rota¬ 
tion curve and an /! 3 -y correlation over the projected length 
of the bar (Bureau & Athanassoula 2005). Instead, even the 
unsymmetrized kinematics clearly show a smooth increase in 
the radial velocity from the southeast side of the galaxy to 
the northwest side, and that is anti-correlated with V, as is 
expected for axisymmetric systems. 

5. MEASURING THE PSF 

The PSF of the NIFS and PPAK observations are important 
inputs into the stellar dynamical models. In order to estimate 
these quantities, we convolve the MGE model presented in 
Section 3 with the sum of two concentric, circular 2D Gaus- 
sians in order to match the collapsed NIES and PPAK data 
cubes. The PSP is parameterized by the dispersion and rela¬ 
tive weight of each Gaussian component. We find dispersions 
of 0"16 and 0."43 with relative weights of 0.61 and 0.39, re¬ 
spectively, for the NIPS PSP, while the PPAK PSP can be 
described with dispersions of 1752 and 5."45 with relative 
weights of 0.82 and 0.18. In addition, the comparison be¬ 
tween the HST image and the collapsed data cubes allows for 
the center of the NIPS and PPAK apertures to be defined. The 
core of the NIPS PSP is larger than expected for AO obser¬ 
vations (e.g., Krajnovic et al. 2009; Seth et al. 2014), which 
is likely the result of using a fairly faint, off-axis tip-tilt star. 
In Section 7.1, we test the effect of our assumed NIPS PSP 
on the inferred black hole mass by instead estimating the PSP 
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TABLE 2 
NIFS Kinematics 


x(") 

(1) 

yi”) 

(2) 

V (km s ' ) 
(3) 

AV (km s * ) 

(4) 

CT (km s ^) 

(5) 

Act (km s *) 
(6) 

h 

(7) 

A/n 

(8) 

h4 

(9) 

A/J4 

(10) 

-0.027 

0.022 

-8.879 

7.817 

396.368 

10.793 

0.017 

0.016 

0.029 

0.017 

-0.027 

-0.028 

-38.385 

6.554 

394.762 

8.844 

0.017 

0.014 

0.033 

0.018 

0.055 

0.022 

37.125 

5.608 

395.394 

7.818 

-0.015 

0.012 

0.031 

0.015 

0.055 

-0.028 

11.045 

6.073 

392.972 

8.368 

-0.001 

0.012 

0.029 

0.014 

-0.077 

-0.003 

-45.778 

6.763 

392.601 

9.207 

0.026 

0.014 

0.028 

0.017 


Note. — Table 2 is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its 
foiTTi and content. Columns (I) and (2) are the x and y Voronoi bin generators, measured relative to the galaxy center. Columns (3) 
- (8) provide the bi-symmetrized NIFS kinematics and en'ors. The position angle is 141.14°, measured counter-clockwise from the 
galaxy’s major axis to x. 
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Fig. 3.— Shown in the top, middle, and bottom panels, respectively, are example PPAK spectra extracted from three spatial locations. Overplotted in red is the 
optimal stellar template convolved with the best-fitting LOSVD, and the gray shaded boxes denote the wavelength regions excluded during the spectral fitting, 
due to the presence of emission lines or sky lines. The model residuals are shown in green, and have been shifted by an arbitrary amount. 


TABLE 3 

PPAK Kinematics 


x(") 

(1) 

y(") 

(2) 

V (kms“‘) 
(3) 

AV (km s ' ) 

(4) 

CT (km s ^) 

(5) 

Act (km s *) 
(6) 

h 

(7) 

A/n 

(8) 

/Z4 

(9) 

A/J4 

(10) 

-0.063 

-0.208 

-28.540 

4.759 

296.999 

6.723 

0.010 

0.013 

0.027 

0.016 

-0.063 

0.792 

18.003 

4.564 

294.131 

6.221 

-0.013 

0.013 

0.025 

0.016 

0.937 

-0.208 

60.969 

4.491 

282.870 

5.858 

-0.030 

0.013 

0.035 

0.017 

-1.063 

-0.208 

-85.825 

4.277 

278.884 

5.984 

0.043 

0.012 

0.036 

0.017 

-0.063 

-1.208 

-60.969 

4.491 

282.870 

5.858 

0.030 

0.013 

0.035 

0.017 


Note. — Table 3 is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its 
form and content. Columns (1) and (2) are the x and y Voronoi bin generators, measured relative to the galaxy center. Columns (3) 
- (8) provide the bi-symmetrized PPAK kinematics and enors. The position angle is 140.71°, measured counter-clockwise from 
the galaxy’s major axis to a;. 
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from NIFS observations we acquired of the tip-tilt star itself. 

6. ORBIT-BASED MODELS 

Black hole masses are often measured from stellar kine¬ 
matics by constructing dynamical models based upon the 
Schwarzschild superposition method (Schwarzschild 1979), 
and here we use the three-integral, triaxial Schwarzschild 
code of van den Bosch et al. (2008). This technique finds 
a self-consistent distribution function from the observables 
without any assumptions about the orbital anisotropy. In the 
model, the black hole, the stars, and dark matter all contribute 
to the galaxy’s gravitational potential. The stellar potential is 
determined by deprojecting the observed surface brightness 
assuming a viewing orientation and a stellar mass-to-light ra¬ 
tio (T) that is constant with radius. We then generate a rep¬ 
resentative orbit library in the potential, and the orbits are 
numerically integrated while keeping track of their intrinsic 
and projected properties. During the modeling, the effects of 
the PSF and aperture binning are taken into account. Finally, 
we assign weights to each orbit such that the superposition 
matches the observed kinematics and total light distribution. 
We calculate many models varying the parameters of inter¬ 
est (Mbh, T, the viewing orientation parameters, and the dark 
matter halo parameters), and the best-fit model is the one with 
the lowest x^. 

6 . 1 . Application to NGC 1271 

For NGC 1271, we adopt a (nearly) oblate axisymmetric 
shape, with an intermediate to long axis ratio of 0.99. We 
assume axisymmetry given that NGC 1271 looks highly flat¬ 
tened and exhibits rapid rotation, and we do not find any evi¬ 
dence for kinematic twists. With this assumption, the inclina¬ 
tion angle (/) is the only viewing orientation parameter needed 
to describe the galaxy’s intrinsic shape. In our final model, we 
adopt i = 83° (where i = 90° corresponds to an edge-on view). 
The MGE description of NGC 1271 presented in Section 3 in¬ 
cludes a couple of flat Gaussian components, which can only 
be deprojected for angles 77 < / < 90°. Therefore, we chose 
to run models at a fixed inclination angle that lies midway be¬ 
tween the extremes. Our choice of inclination angle is further 
supported by a nuclear dust disk that appears highly inclined 
in an HST WFC3 F814W image. While the F814W image 
clearly shows a regular dust lane, the F160W image and the 
/f-band NIFS data cube used in our analysis do not appear to 
be significantly affected by dust. In addition, the NGC 1271 
models include a spherically symmetric dark matter halo fol¬ 
lowing a Navarro-Frenk-White (NFW) form (Navarro et al. 
1996). The parameters describing the NFW halo are the con¬ 
centration index (c) and the fraction of dark matter (/dm), 
where /dm = and Mdm is the halo virial mass and 

M* is the stellar mass. Hence, our models have four free pa¬ 
rameters: Mbh, Th, c, and /dm- 

We began by calculating models on a coarse grid that 
spanned a wide range of values before generating models on 
a smaller, more finely sampled grid focused around the mini¬ 
mum x^- Ultimately, our final model grid contained 21 Mbh 
values, 31 T// values, 8 c values, and 20 /dm values with 
log(MBH/M 0 ) G [9.2,10.2], Th (Tq) G [0.5,2.0], c G [2,16], 
and /dm G [50,1000]. The orbit library samples 29 equipo- 
tential shells with radii between 0"003 to 100", and 8 angular 
and 8 radial values at each energy. We note that triaxial orbit 
families (e.g., box orbits) continue to be included in our or¬ 
bital libraries because we are running a triaxial Schwarzschild 


code in the axisymmetric limit. Moreover, we employ a 
dithering method, in which 125 orbits with adjacent initial 
conditions are bundled together, to establish a smooth dis¬ 
tribution function when constructing Schwarzschild models. 
Thus, the galaxy models are made with 696,000 orbits. The 
models were fit to the observed NIFS and PPAK kinematics, 
where 4 GH moments were measured in a total of 395 bins, 
resulting in 1580 observables. 

7. MODELING RESULTS 

The results of the final model grid are summarized in Fig¬ 
ure 4, which shows the as a function of Mbh, c, and 
/dm after marginalizing over the other three parameters. We 
find best-fit values of Mbh = 3.0 x 10® Mq, Th = 1.40 T 0 , 
c= 16, and /dm = 50 (corresponding to Mdm = 5.Ox 10'^ Mq). 
The comparison between the observed NIFS kinematics and 
model predictions is shown in Figure 5 and the comparison 
between the PPAK data and model predictions is displayed in 
Figure 6 . The model is an excellent match to the observed 
kinematics, and is able to reproduce the sharp increase in the 
velocity dispersion and the slight peak in /i 4 at the nucleus. 
The x^ per degree of freedom (x^) is 0.3. If the unsym¬ 
metrized kinematics are used instead, the best-fit model has 
Xy = 1.3 and an Mbh and Th that are consistent within the 
final uncertainties given below in Section 7.1. 

As can be seen in Figure 4, we are able to place strong 
constraints on Mbh and Th. By searching for the range 
of Mbh, or Th, values that caused the minimum to in¬ 
crease by 1 and 9, we estimate the lu and 3 (j model fitting 
uncertainties. We find Mbh = (3.0^Qg) x 10® Mq (Ict) and 
Mbh = (3.0!?;^) x 10® Mq (3cr), as well as Th = 1-40!|;;J]^ Tq 
(Icr) and Th = 1.40^g [5 Tq (3a). In contrast, the dark halo 
parameters are not well constrained, and the c and /dm val¬ 
ues are highly uncertain. Specifically, the marginalized x^ 
curve for c is unconstrained at the upper end, with c > 5 at 
the 3a level. Although the best-fit value for c is at the maxi¬ 
mum value considered in our final model grid, previous runs 
using coarsely sampled grids showed no sign of convergence 
for large values of c. Similarly, the marginalized y^ curve for 
/dm is essentially flat, and all /dm values sampled by our grid 
(50 < /dm < 1000) are allowed within 3a. The best-fit value 
found for /dm is the minimum value considered in our final 
grid, but we also ran models without a dark halo and found a 
significantly worse fit, such that the y^ increased by 123 rela¬ 
tive to the best-fit model with a dark halo. Thus, a dark halo is 
required to match the observed kinematics, but its properties 
cannot be pinned down with the current datasets and we do 
not present any further details associated with the dark halo 
parameters. 

7.1. Error Budget 

The formal model fitting uncertainties quoted in the previ¬ 
ous section are the statistical errors associated with the dy¬ 
namical models, however systematic effects can have a sig¬ 
nificant impact on the inferred black hole mass and mass-to- 
light ratio. In this section, we evaluate some common sources 
of uncertainty that are not already incorporated into the statis¬ 
tical errors, such as those associated with the assumed form of 
the dark matter halo, the adopted inclination angle, the num¬ 
ber of orbits used in the models, the NIFS PSF model, details 
associated with the extraction of the NIFS kinematics, and the 
symmetrization of the input kinematics. 
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IosCMbh/Mo) 


Th (To) 



Fig. 4.— The plots summarize the results of the stellar dynamical models run with i = 83° and an NFW dark matter halo. The is shown as a function of 
3 ^bh (top left), T/r (top right), c (bottom right), and /dm (bottom left) after marginalizing over the other three parameters. The red point denotes the best-fit 
model, and the dashed line depicts where minimum x^ has increased by 9, which corresponds to the statistical 3(T uncertainties. The uncertainty in Mbh and TT H 
due to systematic effects is significantly larger than that statistical Icr error (see Section 7.1). Thus, we do not plot the statistical Itr confidence level, as it is not 
representative of the range of possible parameter values. The A/bh and T h parameters are well constrained despite the large uncertainties associated with the 
dark halo parameters. 


Dark Matter Halo: Previous work has clearly shown that 
the dark halo could be an important component in the stel¬ 
lar dynamical models due to the degeneracy between the dark 
halo and stellar mass-to-light ratio, which in turn is also de¬ 
generate with the black hole (e.g., Gebhardt & Thomas 2009; 
Schulze & Gebhardt 2011; Rush et al. 2013). If rsphere is spa¬ 
tially well resolved, then the degeneracy between the black 
hole and mass-to-light ratio can be mitigated, and the exclu¬ 
sion of a dark halo will have minimal impact on the inferred 
Mbh- In contrast, if rsphere is not very well resolved and a dark 
halo is not included in the stellar dynamical models, then T 
will be artificially elevated to account for the missing mass 
at large radii. Since T is taken to be constant with radius, a 
smaller Mbh is then required to fit the observed central kine¬ 
matics. 

When generating a grid of models without a dark halo for 
NGC 1271, and assuming i = 83° to match our fiducial model 
presented in Section 7, we measure Mbh = 1.4 x 10^ Mq and 
Th= 1-85 Tq. In other words, the inferredMBH is underesti¬ 
mated by a factor of '^2 when a dark halo is excluded from the 
modeling. We note that when fitting dynamical models with¬ 
out a dark halo to only the small-scale NIFS kinematics, we 


recover very similar parameter values as those found from our 
fiducial model, with Mbh = 2.9 x 10® Mq and T// = 1-45 T 0 . 
Moreover, we tested how the assumed shape of the dark halo 
affects Mbh- Our fiducial model was calculated assuming 
a spherical NFW halo, but another common form is a halo 
with a cored logarithmic potential (Binney & Tremaine 1987; 
Thomas et al. 2005), given by 


Pdm(t) - 


3r^ + r^ 
47rG (r^ + r^)^ 


( 1 ) 


The parameters Vc and r^ are the asymptotic circular velocity 
and radius within which the dark matter density is constant. 
Thus, the halo from a logarithmic potential yields smaller 
densities at small radii compared the the NFW halo. We con¬ 
structed models with a dark halo from a cored logarithmic po¬ 
tential, sampling 9.2 < log(MBH/MQ) < 10.2,0.5 < T// < 2.0 
Tq, 100 <Vc < 700 km s“', and 1 < < 32 kpc. We recov¬ 

ered similar results to the models with an NFW halo, namely 
that Mbh = 3.2 x 10® Mq and Th = 1-35 Tq, corresponding 
to a 7% increase in Mbh and a 4% decrease in T// compared 
to the fiducial model. Although the black hole mass is sen- 
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Fig. 5.— The bi-symmetrized NIFS kinematics for NGC 1271 (top) are compared to the best-fit model predictions (bottom), where Mbh = 3.0 x 10^ Mq and 
T/f = 1.40 T 0 . The same scaling, shown by the color bar on the right with the minimum and maximum values given at the top of the maps, is used to plot the 
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Fig. 6 .— The bi-symmettized PPAK kinematics (top), which extend out to about 4 Re, are shown along with the predictions from the best-fit stellar dynamical 
model (bottom); see Figure 5 for description. Kinematic measurements are missing from x~-15”, y~10” and from x~20”, y~0” due to the presence of 
foreground objects, which were masked before extracting the kinematics. The southeast side of the galaxy has blue-shifted radial velocities. 


sitive to the inclusion of a dark halo in the stellar dynamical 
models, the form of the halo has a small impact on Mbh- For 
this reason, being unable to constrain the dark halo parame¬ 
ters is not a concern for the purposes of this paper, as long 
as reasonable halos are sampled over when constructing the 
orbit-based models. 

Inclination Angle: All of the models presented in this pa¬ 
per assume an axisymmetric shape with an inclination angle 
of i = 83°. However, we also ran a grid of models that sam¬ 
pled 13 inclination angles from 77 to 89°. This corresponds 
to the range of angles for which the MGE model in Section 
3 can be deprojected. It is computationally expensive to cal¬ 
culate a model grid that samples over Mbh, "f/r, i, c, and /dm 


simultaneously, so we instead varied the first three parame¬ 
ters while sampling over five NEW halos. The five dark halos 
were those with the lowest from the model grid at the be¬ 
ginning of Section 7, and the halos span a range of c and /dm 
values (8 < c < 16 and 50 < /dm < 500). Erom this test, we 
determine that Mbh = 3.3 x 10^ Mq, which is within 10% of 
the best-fit value in Section 7, and T// = 1.40 T 0 , which is the 
same as the best-fit value in Section 7. Moreover, the best-fit 
inclination angle was / = 87°, however the angle was not well 
constrained. All angles between 79° and 89° were allowed 
within the 3cr statistical uncertainties. Such behavior is not 
surprising, and other stellar dynamical work have also found 
it difficult to infer the inclination angle from 2D line-of-sight 
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kinematics (Krajnovic et al. 2005; van den Bosch & van de 
Ven 2009). 

Number of Orbits: The fiducial model presented in Sec¬ 
tion 7 was calculated using orbits that covered 29 equipoten- 
tial shells with 8 angular and 8 radial values at each energy. 
When accounting for the orbital dithering, this translates into 
a total of 696,000 orbits. We also tested the effect on Mbh 
and Th when the number of orbits is about doubled, such that 
37 equipotential shells with 10 angular and 10 radial values 
at each energy are used. Again 125 orbits with adjacent start¬ 
ing positions were bundled together, resulting in 1,387,500 
orbits. Due to a large increase in computational time for a 
single model, we constructed a model grid that samples over 
Mbh, ^h, and the top five NFW halos from Section 7. We 
found no change in the best-fit Mbh or T values compared to 
the fiducial model. 

NIFS PSF: The NIFS PSF was measured by comparing the 
MGE model of the HST image to the collapsed NIFS data 
cube. While this approach is commonly used in black hole 
mass measurement work (e.g., Krajnovic et al. 2009; Seth 
et al. 2010; Walsh et al. 2012), estimation of the PSF from 
AO observations is notoriously difficult due to the constantly 
changing quality of the AO correction and the combination of 
data cubes from multiple nights. Therefore, we also estimated 
the PSF in a different manner in order to assess how strongly 
the adopted NIFS PSF affects Mbh. Utilizing the NIFS obser¬ 
vations of the tip-tilt star, we fit the sum of four concentric, 
circular 2D Gaussians to the collapsed NIFS data cube using 
Galfit. We found that the PSF is best described by Gaussians 
with dispersions of 0"04, 0"08, 0."21, and 0"43 with relative 
weights of 0.08, 0.36, 0.26, and 0.30. A four-component PSF 
model provided a significantly better fit to the collapsed data 
cube than a simpler two-component model, and the residuals 
between the four Gaussian model and the data have a standard 
deviation of just 9% out to a radius of 1". The PSF is typical 
of what one would expect from the Gemini AO system (e.g., 
Gebhardt et al. 2011; Onken et al. 2014), and the quality of the 
AO correction is better than that implied by the PSF adopted 
in the fiducial model. This is not completely surprising as 
the NIFS observations of the star used the star itself for tip- 
tilt corrections, whereas the observations of the galaxy were 
made off-axis, using the star for guiding. Nonetheless, cal¬ 
culating stellar dynamical models using this better PSF and 
comparing to the results using the poorer PSF should cover 
the range of possible black hole masses due to the uncertainty 
in the NIFS PSF. We calculated models using the new four- 
Gaussian PSF and further assumed that the center of the NIFS 
spaxel with the largest flux coincides with the galaxy nucleus. 
We varied Mbh and T^, while sampling over the top five 
NFW halos from Section 7, and found Mbh = 2.7 x 10® Mq 
and Th= 1 45 Tq. Therefore, the black hole mass and mass- 
to-light ratio change by 10% and 4% compared to the fiducial 
values in Section 7. 

Measuring the NIFS Kinematics: We measured the NIFS 
kinematics with pPXF using a second degree additive Leg¬ 
endre polynomial and a second degree multiplicative polyno¬ 
mial to correct for shape differences between the LOSVD- 
broadened optimal stellar template and the observed galaxy 
spectrum. We selected this continuum correction because it 
was one of the simplest models that still provided a good fit 
to the data, and produced kinematics that were in good agree¬ 
ment with those measured using combinations of degree 0-3 
additive/multiplicative polynomials, with the exception of the 


lowest order polynomials. Since the continuum has been pre¬ 
viously divided out of the the spectral templates in the NIFS 
library, but not removed from the observed galaxy spectra, 
using the lowest order polynomials (combinations of degree 
0-1 additive/multiplicative polynomials) produced visibly 
poor spectral fits. 

We also constructed a new library containing stellar tem¬ 
plates that have not been continuum-corrected and thus have a 
very similar shape as the galaxy spectra. We retrieved X-band 
NIFS observations from the Gemini archive of stars that are 
part of the NIFS Spectral Template Library and two stars ob¬ 
served under program GN-2010A-Q-112. These twelve stars 
are KO - M5 giants, a K5 supergiant, and an MO supergiant. 
We reduced the observations following the main procedure 
outlined in Section 2.2, with the additional steps of extracting 
a one-dimensional spectrum, rebinning to a common wave¬ 
length range and sampling, and shifting the stars to rest. With 
this new template library, we are able to obtain good fits to the 
galaxy spectra using low-order polynomials with pPXL. 

In order to examine possible effects on Mbh and Th due 
to uncertainties associated with the choice of the pPXL poly¬ 
nomial degree, we fit dynamical models to the NIPS kine¬ 
matics extracted using the new stellar template library and an 
additive constant, along with the PPAK kinematics presented 
in Section 4.2. We sampled over Mbh, T h, and the top five 
NEW halos from Section 7, finding Mbh = 2.4 x 10® Mq and 
Uh = 1.45 Tq. This corresponds to a change of 20% and 
4% from the best-fit black hole mass and mass-to-ratio from 
the fiducial model. These results are likely representative of 
a maximum change in best-fit parameter values, as using this 
particular continuum correction with the new stellar template 
library resulted in the largest number of bins with kinematics 
inconsistent at the Icr level (all bins were consistent at the 2a 
level) compared to the NIPS kinematics from Section 4.1. 

Symmetrizing the Kinematics: The NIPS and PPAK kine¬ 
matics were bi-symmetrized prior to using them as inputs into 
the stellar dynamical modeling code. We used the method 
outlined in van den Bosch & de Zeeuw (2010), which aver¬ 
ages the measurements of a GH moment in a four-fold sym¬ 
metric manner around the minor and major axes. During the 
averaging of the measurements for a single GH moment, the 
kinematic error of the bin and the fraction of flux a spaxel con¬ 
tributes to that bin are taken into account. Such modifications 
to the input kinematics and their errors are a common way in 
which to reduce observational noise, and often kinematics are 
bi-symmetrized in order to obtain reasonable results from an 
axisymmetric modeling code and are point-symmterized for 
use with a triaxial modeling code. However, symmterization 
routines are never perfect, as discussed in van den Bosch & 
de Zeeuw (2010). We therefore tested how symmetrization 
affects Mbh and T// by running additional models with kine¬ 
matics that were point-symmetrized. We find that the black 
hole mass increases to 3.6 x 10® Mq, or by 20% of the fidu¬ 
cial value, and the //-band mass-to-light ratio decreases to 
1.35 Tq, or 4% of the fiducial value from Section 7. 

Summary: We derive the final range of black hole masses 
and mass-to-light ratios for NGC 1271 by adding in quadra¬ 
ture the formal model fitting Icr uncertainty from Section 7 
and the additional sources of systematic uncertainty above. 
Ultimately, we determine that Mbh = (3.0);[ j) x 10® Mq and 
Th = 1.40^g[j Tq. The dominant source of systematics 
for NGC 1271 are those associated with the continuum- 
correction model used to extract the NIPS kinematics and the 
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symmetrization of the kinematics. Both affect Mbh and T// at 
the 20% and 4% level, respectively. 

7.2. Other Considerations 

In addition to examining possible sources of systematic un¬ 
certainty and incorporating the effects into the final error bud¬ 
get, we ran other tests to assess the robustness of the NGC 
1271 Mbh measurement. We describe these tests below. 

Nuclear Dust Disk: A dust disk is present at the center of 
NGC 1271, and is visible in the F814W WFC3 image. How¬ 
ever, dust doesn’t appear to be significant in the near-infrared 
imaging and in the NIFS data cube. We tested constructing 
a new MGE of the F160W image after excluding the dust 
disk by using the F814W image as guide for generating the 
mask. We ran dynamical models using the modified MGE, 
following the procedure in Section 6.1 and sampling over the 
top five NEW halos from Section 7. We found best-fit val¬ 
ues of Mbh = 3.2 x 10® Mq and Th = 135 Tq, which is well 
within the final uncertainties adopted for NGC 1271 discussed 
in Section 7.1. 

Stellar Mass-to-Light Ratio Variation: Our dynamical mod¬ 
els assume that Th remains constant with radius. In order to 
determine whether there is an obvious change in stellar popu¬ 
lation, we generated an MGE of the E814W image following 
the methods described in Section 3. During the fit, we masked 
out the nuclear dust disk and foreground objects in the E814W 
image, and we accounted for the PSE using a bright, isolated 
star in the image. Erom the MGE models of the E814W and 
E160W images, we don’t see evidence for color gradients, 
finding that the color changes by at most 0.16 mag from 0"2 
to 24". 

Although NGC 1271 exhibits a fairly uniform color, we 
further examined dynamical models that are fit to only the 
NIES kinematics. Given the limited radial extent of the NIES 
kinematics, which extend out to a radius of ~ 1 ", or ^390 
pc, systematics associated with mass-to-light ratio gradients 
(and dark matter halos) are mitigated. When fitting to only 
the NIES kinematics, we recover consistent results to those 
presented in 7.1, where Mbh = (3.5^^o) x 10® Mq and Th = 
1.30;!^ q 7 Tq (1(T uncertainties). We show contours of as 
a function of black hole mass and mass-to-light ratio for this 
NIES-only model grid in Eigure 7. 

PPAK Kinematics: When measuring the stellar kinematics 
from the PPAK data, we masked out spectral regions con¬ 
taining possible emission lines and sky residuals, as can be 
seen in Eigure 3. In order to verify that our choice of a spec¬ 
tral mask does not bias the kinematic measurements and in¬ 
ferred black hole mass, we decreased the number and width 
of the excluded wavelength regions. We ran dynamical mod¬ 
els using the modified PPAK kinematics while sampling over 
the top five dark matter halos in Section 7. We measured 
Mbh = 3.2 x 10® Mq and Th = 1.35 Tq, which is within the 
final uncertainties given for NGC 1271 in Section 7.1. 

Also, the observed velocity dispersion of the PPAK kine¬ 
matics presented in Section 4.2 drops below the instrumen¬ 
tal resolution in many of the spatial bins located > 10 " away 
from the nucleus. While care was taken to homogenize the 
line spread function to a common value such that there was 
no variation with wavelength or fiber position, measuring dis¬ 
persions well below the instrumental resolution is a difficult 
task. We therefore also tested the effect on Mbh and Th when 
excluding the spatial bins in which the dispersion was below 
150 km s“'. When calculating dynamical models that fit to 
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Fig. 7. — Contours of shown as a function of black hole mass and 

//-band mass-to-light ratio for the case when dynamical models are fit to just 
the NIFS kinematics. At each gray point a model was calculated, and the red 
square denotes the best-fit model. The red contour and two black contours 
signify where has increased by 1, 4, and 9 from the minimum. Thus, the 
vertical lines show the Icj uncertainties on Mbh and the horizontal lines give 
the Icj uncertainties for T// when marginalizing over the other parameters. 

the adjusted PPAK kinematics and that sample over the top 
five dark matter halos in Section 7, we find no change from 
the best-fit values presented in Section 7. 

8 . DISCUSSION 

NGC 1271 harbors a black hole with Mbh = (3.0!i;?) x 
10® Mq and has a stellar mass-to-light ratio of Th = 
1.40tl? Tq. We note that the final uncertainty on the black 
hole mass we use is comparable to the formal 3a statistical 
uncertainty. Some (Cappellari et al. 2009; Krajnovic et al. 
2009; Emsellem 2013) have suggested 3cr statistical errors 
should be used in place of Icr errors as a conservative way 
in which to account for the effect of unknown systematics on 
Mbh- With a black hole mass of 3.0 x 10® Mq and adopting 
276 km s“' for the bulge stellar velocity dispersion (see Sec¬ 
tion 8.2), Tsphere = 0"44. Thus, the NIPS observations have 
resolved the black hole sphere of influence. Below we dis¬ 
cuss the galaxy’s orbital structure and place the galaxy on the 
MBH-host galaxy relations. 

8 .1. Orbital Structure 

In addition to determining the mass of the black hole in 
NGC 1271, the Schwarzschild models provide information 
about the galaxy’s orbital structure. Using our best-fit model 
in Section 7, we show the ratio Urja, as a function of radius 
in Eigure 8 . The tangential velocity dispersion is defined as 
af = {a^ + ag)/!, and (r,9,(j>) are the usual spherical coor¬ 
dinates. We find that NGC 1271 is roughly isotropic at all 
radii covered by our kinematic measurements, deviating by at 
most 30% from arjat = 1, but we observe a trend in which 
Of!a, declines at radii outside the black hole sphere of in¬ 
fluence. As expected, short-axis tube orbits dominate in this 
oblate system, making up more than 85% of the orbits at all 
radii. Long-axis tube orbits, which are important for triaxial 
and prolate systems, are negligible, while the fraction of box 
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Fig. 8 .— NGC 127rs orbital structure, as inferred from the best-fit dy¬ 
namical model, is shown. The anisotropy (top) and orbit type (bottom) are 
displayed with radius over the range covered by the NIFS and PPAK kine¬ 
matic measurements. The horizontal dashed gray line denotes the isotropic 
case and the vertical dot-dashed gray line shows the black hole sphere of 
influence. NGC 1271 is roughly isotropic at all radii and is dominated by 
short-axis tube orbits as is expected for oblate systems. 


orbits increases at small radii but still make up only 15% of 
the orbits near the nucleus. 

Furthermore, we use our best-fit stellar dynamical model to 
examine the mass distribution as a function of average radiu^, 
F, and spin, A^, of the orbits, where X^=f^x (jja)- Here, 
is the average angular momentum along the z-direction and a 
is the average second moment of the orbit. NGC 1271 shows 
several dynamical components, as can be seen in Figure 9, 
including a clear non-rotating bulge (with -0.2 < A^ < 0.2), 
a highly co-rotating component (with Aj_^0.5), and a maxi¬ 
mally co-rotating component (with 0.8 < Aj < 1.0). The bulge 
component accounts for 12% of the mass, whereas the rotat¬ 
ing components total 75% of the mass within the radial ex¬ 
tent of the kinematic measurements. Qualitatively these com¬ 
ponents agree with the classification of NGC 1271 as a fast- 
rotating SO galaxy with a classical bulge. 

8.2. Black Hole - Host Galaxy Relations 

Although the dynamical decomposition from the best-fit 
stellar dynamical model presented above can be used to place 
NGC 1271 on the Mbh - bulge relationships, we follow the 
more common approach of carrying out a photometric de¬ 
composition to determine the galaxy’s bulge luminosity and 
bulge effective radius. Using Galfit, we find a single Sersic 
component fit is an insufficient description of the galaxy, with 
the percent difference between the model and data reaching 
as high as 60%. The fit is significantly improved with the 
addition of one or two other Sersic components, and in the 
later case the percent difference between the model and data 
is under 15%. In Table 4, we present the best-fit parameters 
of Galfit models with one, two, and three Sersic components, 
as well as the F160W luminosity for each component. The 
three-component model provides the best match to the HST 
image, but it is difficult to unambiguously identify a “bulge” 
component because the components all have rather low Ser¬ 
sic indices. Therefore, we conservatively assume that the in¬ 
nermost component of the three-component model provides 
a lower limit on the bulge luminosity and effective radius, 
while the innermost component of the two-component model 
gives an upper limit. This yields a /f-band bulge luminosity 
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Fig. 9.— The mass distribution is plotted as a function of average spin and 
radius of the orbits for the region covered by our kinematics measurements. 
The dynamical decomposition utilizes the orbital weights from our best-fit 
stellar dynamical model and shows distinct non-rotating bulge (- 0.2 < < 

0 . 2 ) and rotating (A^ > 0 . 2 ) disk components. 


TABLE4 
Galfit Models 


Component 

Mh 

Lh (Lq) 

Re (") 

n 

g' 

(1) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

1 

10.72 

7.7xl0‘° 

5.56 

4.78 

0.41 

1 

10.95 

6 . 2 xl 0 ‘° 

5.23 

6.54 

0.54 

2 

12.47 

1.5x10“* 

6.63 

0.93 

0.20 

1 

12.39 

1 . 6 x 10 “* 

0.61 

2.12 

0.68 

2 

11.80 

2 . 8 x 10 “* 

5.24 

1.17 

0.27 

3 

11.90 

2 . 6 x 10 “* 

10.46 

1.30 

0.62 


Note. — Column (I) .shows the Sersic component number, column (2) provides the 
FI60W apparent magnitude in the Vega system, not yet corrected for galactic extinction, 
column (3) gives the luminosity after a correction of 0.085 for galactic extinction and 
assuming an absolute solar magnitude of 3.33, column (4) is the effective radius, column 
(5) gives the Sersic index, and column (6) lists the projected axis ratio. 

of (1.9-7.2) X 10'° Lq when correcting for galactic extinc¬ 
tion using the Schlafly & Finkbeiner (2011) WFC3 F160W 
value of 0.085, assuming an H-K color of 0.2 (Vazdekis 
et al. 1996), and a /T-band solar absolute magnitude of 3.29. 
Alternatively, the bulge has a mass of (2.3-8.7) x 10'** Mq 
when applying the best-fit mass-to-light ratio from our dy¬ 
namical models to the luminosities in Table 4. The corre¬ 
sponding bulge effective radius ranges between 0."6 and 5"2 
(or 0.2-2.0 kpc). 

In order to determine the bulge stellar velocity dispersion 
for NGC 1271, we use the best-fit stellar dynamical model 
from Section 7 and predict the luminosity-weighted second 
moment within a circular aperture whose radius equals the 
galaxy’s bulge effective radius, following the approach used 
by van den Bosch et al. (2012). Due to the uncertainty in the 
bulge effective radius for NGC 1271, we measure the effec¬ 
tive stellar velocity dispersion (cr^^bui) for three different bulge 
effective radii corresponding to the largest estimate from 
the Galfit decomposition, the smallest estimate of R^, and the 
average of the two. Additionally, some previous black hole 
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studies have chosen to exclude data within rsphere when de¬ 
termining (Te^bui because the stellar kinematics are under the 
direct influence of the black hole in this region (e.g., Geb- 
hardt et al. 2011; McConnell & Ma 2013). When exclud¬ 
ing the region within rsphere, we And effective stellar velocity 
dispersions of (Te_bui = 272 km s“', (Te.bui = 276 km s“', and 
o'e,bui = 349 km s“', whereas when the region within rsphere is 
included we measure CTg^bui = 285 km s“', ag^hui = 294 km s“', 
and (Te^bui = 358 km s“', for bulge effective radii of 5."2, 2."9, 
and 0."6, respectively. As a comparison, the HET Massive 
Galaxy Survey reports a central velocity dispersion of 317 km 
s“' (van den Bosch et al. 2015), which is the observed stellar 
velocity dispersion within a 3 "5 aperture from the major axis 
long-slit data. 

Figure 10 shows the location of NGC 1271 on the most 
recent versions of the Mbh - o'* and Mbh - Cbui relations by 
Kormendy & Ho (2013). For the purposes of placing NGC 
1271 on the Mbh - o'* correlation, we use ct* = 276 km s“' 
with uncertainties that include the o'e^bui measurements made 
for bulge effective radii of 5."2 and 0"6 when excluding data 
within rsphere- When placing NGC 1271 on the Mbh - Cbui 
relation, we set the faint end of the bulge luminosity error 
bar assuming the luminosity of the innermost component of 
the three-component Sersic fit to the HST image and the high 
end of the error bar assuming the luminosity of the innermost 
component of the two-component Sersic fit. We adopt a K- 
band bulge luminosity of 4.6 x 10 '° Lo, which is the midpoint 
of the range of possible bulge luminosities. 

We find that NGC 1271 consistent with the Mbh - o'* rela¬ 
tion, but is an order of magnitude above the black hole mass 
prediction from the Mbh - Cbui correlation. In order to demon¬ 
strate that the black hole in NGC 1271 must be larger than 
that expected from MBH-Cbui, in Figure 11 we present the 
NIFS observations of the velocity dispersion and h 4 along 
with the predictions from the best-fitting model with a black 
hole mass of 3.0 x 10° Mq and a 4.7 x 10* Mq black hole. 
The 4.7 X 10* Mq black hole is expected from Mbh - Cbui 
when conservatively using the galaxy’s total /f-band luminos¬ 
ity of 8.9 X 10'° Lq, which is derived from the single Ser¬ 
sic fit to the HST image, after correcting for galactic extinc¬ 
tion and assuming a H-K = 0.2. Clear differences between 
the kinematic predictions and the observations can be seen by 
eye. The best-fit model with Mbh = 3.0 x 10° Mq is able to 
nicely reproduce the sharp rise in the velocity dispersion and 
the slight peak in h 4 at the nucleus, while the less massive 
black hole predicted from Mbh - Tbui fails to do so. 

NGC 1271 has an apparent ellipticity of e = 0.6 
and a specific stellar angular momentum of Ar = 
{R\V\)/{RVV^ + a^) = 0.5 within one effective radius 
based on the PPAK data. Here, R, V, and cr are the radius, 
velocity, and velocity dispersion and the brackets denote a 
luminosity weighted average (Emsellem et al. 2007). Using 
the dividing line between slow and fast rotators from the 
ATLAS*'^ survey, such that fast rotators have Ar > 0.31 x 
(Emsellem et al. 2011), NGC 1271 falls well within in this 
fast rotator regime. 

NGC 1271 appears similar to the other compact galaxies 
NGC 1277, NGC 1332, NGC 4342, NGC 4486B, and M60- 
UCDl. All six of these early-type galaxies have small sizes, 
are rotating, show large stellar velocity dispersions for their 
luminosities, and have black holes that are too massive for 
their host galaxy’s bulge luminosity. The black holes, how¬ 


ever, are consistent with Mbh - o'* given the intrinsic scatter of 
the relation. In the case of M60-UCD1, Seth et al. (2014) sug¬ 
gest that the ultracompact dwarf galaxy (UCD) was once the 
nucleus of a larger galaxy that has since been tidally stripped 
by the giant elliptical M60, whose center lies at a projected 
distance of just 6.6 kpc away from the UCD. While tidal 
stripping is a natural explanation for the presence of an over- 
massive black hole, in the case of NGC 1271, we do not see 
signs of active stripping in the HST image. The isophotes ap¬ 
pear extremely regular, and no massive galaxies immediately 
neighbor NGC 1271 like in the case of M60-UCD1. NGC 
1271 is ~270 kpc in projection from the BCG of Perseus. 
Further evidence could come from counting the number of 
globular clusters, as the galaxy would be stripped of its glob¬ 
ular clusters first. While NGC 1271 appears not to have been 
stripped with our current data, we cannot rule out an event in 
the distant past. 

Interestingly, the behavior of the compact, high-dispersion 
galaxies being consistent with Mbh - cr* but being large pos¬ 
itive outliers on Mbh - Tbui could be in conflict with recent 
observations of BCGs, which instead may hint that black 
hole mass becomes independent of cr* at high black hole 
mass while the Mbh - Tbui correlation remains unchanged 
at large luminosities (McConnell & Ma 2013; Kormendy & 
Ho 2013). Clearly, more compact, high-dispersion galaxies 
and BCGs/giant ellipticals need to be examined. There could 
be systematic differences in the scaling relations between the 
two types of galaxies, thereby imply different mechanisms 
for black hole growth. Since the compact, high-dispersion 
galaxies like NGC 1271 look similar to the quiescent z ~ 2 
red nuggets, they could be relics that somehow avoided the 
same fate that ultimately produced the giant ellipticals ob¬ 
served today. Perhaps the compact, high-dispersion galaxies 
are left over from an era when the local black hole scaling 
relations did not apply and galaxies instead contained over- 
massive black holes. 

9. CONCLUSION 

To summarize, we obtained AO-assisted Gemini NIFS ob¬ 
servations of NGC 1271 to map out the stellar kinematics on 
scales comparable to the black hole sphere of influence, and 
large-scale IFU data with PPAK, which are useful for con¬ 
straining the galaxy’s stellar mass-to-light ratio and orbital 
distribution. Using an HST WFC3 //-band image along with 
the spectral information, we constructed orbit-based stellar 
dynamical models. We measure Mbh = (3.0^}'°) x 10° Mq 
and Th = 1.40;!;q}[ Tq. The quoted errors combine the Icr 
model fitting uncertainties with some common sources of sys¬ 
tematic uncertainty that affect stellar dynamical models. The 
black hole in NGC 1271 is at the upper end of the black hole 
mass distribution (Mbh > 1 x 10° Mq). Yet, this compact, 
rapidly rotating galaxy, with a high stellar velocity disper¬ 
sion for its luminosity is very different from the giant ellip¬ 
tical galaxies and BCGs that are expected to harbor the most 
massive black holes in the Universe. Such host galaxy envi¬ 
ronments have yet to be widely explored on the Mbh '- host 
galaxy relations. With our mass measurement, we find that 
the black hole is too large for the galaxy’s //-band bulge lu¬ 
minosity of (4. 6 !|°)xlO"'Lo , falling an order of magnitude 
above the expectation from the Mbh - Tbui correlation, but 
the black hole mass is consistent with expectations from the 
Mbh '- cr* relationship assuming cr* = 276^^ km s“'. This 
behavior has also been observed in the few other compact 
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Fig. 10.— NGC 1271 (red filled square) is shown on the black hole - host galaxy relations. The black hole/galaxy measurements (black and gray filled circles) 
and the fitted relations (dot dashed lines) are taken from Kormendy & Ho (2013). The compact galaxies that have existing dynamical Mbh measurements are 
denoted with the red asterisks. These galaxies are generally consistent with Mbh - <r* but are positive outliers from Mbh - f-bul- Kormendy & Ho (2013) did not 
include the measurements shown in gray and light red when fitting the black hole scaling relations. 
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Fig. 11.— The observed velocity dispersion (top) and /14 (bottom) measured from the NIFS data (left) is compared to predictions from the best-fit model with 
^BH = 3.0 X 10^ Mq (middle) and a model with a 4.7 x 10^ Mq black hole (right), which is the mass predicted from Mbh - ^bui when conservatively adopting 
the total galaxy luminosity. When generating the g and predictions for a 4.7 x 10^ Mq black hole, we sample over Th and the top five NFW dark halos from 
the model grid in Section 7, such that the combination of parameters produces a model with the lowest for a black hole mass of 4.7 x 10^ Mq. The data and 
model maps are plotted on the same scale, with the ranges given by the color bar to the right and the minimum and maximum values printed at the side of the 
maps. The best-fit model is able to reproduce the shaip rise in the velocity dispersion and slight peak in h 4 at the center, while the smaller black hole expected 
from the Mbh - f-bui correlation is unable to match either of these features. 
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galaxies that have dynamical black hole mass measurements 
to date. Carrying out more black hole mass measurements 
in similar galaxies using high spatial resolution observations 
from HST and AO is necessary in order to determine if there 
are systematic differences in the black hole scaling relations 
between the large ellipticals/BCGs and these compact, high- 
dispersion galaxies. The compact, high-dispersion galaxies 
could be remnants of the z ~ 2 red nugets that for some rea¬ 
son did not evolve into the largest ellipticals observed today, 
and instead reflect a time when black holes were too large 
for their bulges. More broadly, additional black hole mass 
measurements are needed in order to enlarge and better fill in 
undersampled regions of galaxy parameter space. Targeting 
such a large and carefully selected sample with high spatial 
resolution facilities is a natural step toward gaining a more 
complete census of local black holes and a better understand¬ 
ing of the role that black holes play in galaxy evolution. 
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